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ABSTRACT 



o 



Aims. The nearby, bright, almost completely unreddened Type la supernova 2011fe in MlOl provides a unique opportunity to test 
both the precision and the accuracy of the extragalactic distances derived from SNe la light curve fitters. 
, Methods. We apply the current, public versions of the independent light curve fitting codes MLCS2k2 and SAUr2 to compute the 

• ■ distance modulus of SN 201 Ife from high-precision, multi-color (BVRI) light curves. 

■ Results. The results from the two fitting codes confirm that 201 Ife is a "normal" (not peculiar) and only slightly reddened SN la. New 
P"! unreddened distance moduli are derived as 29.21 ± 0.07 mag (D ~ 6.95 ± 0.23 Mpc, MLCS2k2), and 29.05 ±0.07 mag (6.46 ± 0.21 
Q '. Mpc). 

5_j ■ Conclusions. Despite the very good fitting quality achieved with both light curve fitters, the resulting distance moduli are inconsistent 

■ by 2(T. Both are marginally consistent (at ~ Icr) with the HST Key Froject distance modulus for MlOl. The SAUr2 distance is 
I in good agreement with the recently revised Cepheid- and TRGB-distance to MlOl by Shappee & Stanek. Averaging all SN- and 

Cepheid-based estimates, the absolute distance to MlOl is ~ 6.6 ± 0.5 Mpc. 

, Key words, supernovae: individual (SN 201 Ife) - galaxies: individual (MlOl) 

>: 
o, 

■ 1. Introduction as well as other details like reddening and intrinsic color for 
£2 ' the calibrating objects. This is the major source of the relatively 

. Supernovae (SNe) la are extensively used for deriving ex- s^all, but still existing random and systematic errors that limit 

1^ . tragalactic distances, because in the last two decades they the precision and accuracy of the derived distances (see e.g. 

O 1 turned^out to b^recise^and reliable distance indicators (see iMandel. Naravan & Kirshneii lloiil for further discussion and 

p\| , e.g. iMatheson et al.L 120121 and references therein). Although references) 
T— I • they are not standard candles in the optical (contrary to the 

; widespread statements that they are, which appear frequendy SN 2011fe (aka PTFl Iklv. iNugent et al.L 1201 ih is an excel- 
even in the most recent papers), their light curve (LC) shape cor- lent object in this respect, because this bright (nip^k ~ 10 mag), 
relates with their peak absolute magnitude, making them stan- neai- by (D ~ 6.5 Mpc) SN la was discovered hours after explo- 
5^ ; dardizable (or calibratable) objects. The main advantage of ap- sion (iNugent et all 1201 iLlBloom et al.Ll2012h and suffered from 
plying SNe la for distance measurement is that the method is very litde interstellar reddening {Ay ~ 0.04 mag, 'Nugen t et all 
essentially photometric and does not need spectroscopy, save for 201 1; Patat et al., 201 1). The high apparent brightness allowed 
typing the SN 'is a la. us to obtain accurate, high signal-t o-noise photometry, wh ile the 
The LC shape is c onnected with the peak brightness via very low redshift (z = 0.000804. Ide Vaucouleurs et al.L |199"T1) 
the "Phillips-relation" (iPhillipsl Il993h . which states that SNe of its host galaxy, MlOl (NGC 5457), eliminates the neces- 
la that decline more slowly after maximum have intrinsi- sity of K-corrections for the photometry that otherwise would 
cally brighter peak magnitude, and vice versa. Although there be a major source of systematic errors plaguing the distance 
are attempts t o explain this pheno r nenon based on theoret- determination dHsiao et al.L 120071: iHsiaol 12009 ). The low in- 
ical grounds (iHoeflich et akl Il996t iPinto & Eastmanl |200U terstellar reddening is also a very fortunate circumstance, be- 
iKasen & Wooslevi l2007h . this "peak magnitude - LC proper- cause all complications regarding the handling of the effect 
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ties" relation is still mostly empirical at present. Therefore, of interstellar dust (galactic vs. non-standard reddening, disen- 
the whole procedure of getting distances from the photome- tangling reddening and intrinsic color variation, etc.) are ex- 
try of SNe la relies on empirical calibrations of the SNe la pected to be minimal. Also, the host galaxy, MlOl, has many re- 
peak brightnesses, which need accurate, independent distances cently published distance estimates by various methods includ- 
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Table 1. BVRI magnitudes of SN 201 Ife from the Konkoly 
Observatory, Hungary 



JD° 


B (mag) 


V (mag) 


R (mag) 


I (mag) 


799.3 


14.394 


(.05) 


14.079 (.05) 


13.982 


(.05) 


13.985 


(.05) 


800.3 


13.688 


(.05) 


13.277 (.05) 


13.230 


(.05) 


13.167 


(.05) 


801.3 


12.880 


(.04) 


12.603 (.04) 


12.580 


(.04) 


12.514 


(.04) 


802.3 


12.203 


(.02) 


12.176 (.02) 


12.069 


(.02) 


12.039 


(.02) 


803.3 


11.738 


(.02) 


11.761 (.02) 


11.656 


(.02) 


11.607 


(.02) 


804.3 


11.327 


(.02) 


11.418 (.02) 


11.264 


(.02) 


11.246 


(.02) 


805.3 


11.054 


(.04) 


11.051 (.04) 


10.958 


(.04) 


10.953 


(.04) 


807.3 


10.480 


(.05) 


10.595 (.05) 


10.509 


(.05) 


10.519 


(.05) 


808.3 


10.299 


(.03) 


10.479 (.03) 


10.372 


(.03) 


10.412 


(.03) 


809.3 


10.141 


(.04) 


10.317 (.04) 


10.276 


(.04) 


10.328 


(.04) 


811.3 


9.955 ( 


.14) 


10.179 (.14) 


10.042 


(.14) 


10.219 


(.14) 


815.3 


10.027 


(.06) 


10.094 (.06) 


10.108 


(.06) 


10.309 


(.06) 


816.3 


10.005 


(.07) 


10.113 (.07) 


10.037 


(.07) 


10.383 


(.07) 


817.2 


10.078 


(.08) 


10.081 (.08) 


10.080 


(.08) 


10.370 


(.08) 


818.2 


10.063 


(.05) 


10.049 (.05) 


10.042 


(.05) 


10.448 


(.05) 


819.2 


10.136 


(.07) 


10.078 (.07) 


10.104 


(.07) 


10.525 


(.07) 


820.3 


10.096 


(.06) 


10.020 (.06) 


10.091 


(.06) 


10.485 


(.06) 


821.3 


10.211 


(.05) 


9.995 (.05) 


10.080 


(.05) 


10.493 


(.05) 


822.3 


10.138 


(.07) 


10.058 (.07) 


10.162 


(.07) 


10.562 


(.07) 


826.3 


10.575 


(.06) 


10.330 (.06) 


10.487 


(.06) 


10.879 


(.06) 


828.3 


10.799 


(.02) 


10.505 (.02) 


10.659 


(.02) 


10.921 


(.02) 


829.3 


10.956 


(.02) 


10.566 (.02) 


10.730 


(.02) 


10.918 


(.02) 


830.3 


11.110 


(.05) 


10.630 (.05) 


10.735 


(.05) 


10.934 


(.05) 


831.3 


11.181 


(.03) 


10.654 (.03) 


10.717 


(.03) 


10.873 


(.03) 


832.3 


11.263 


(.02) 


10.744 (.02) 


10.774 


(.02) 


10.856 


(.02) 


835.2 


11.689 


(.02) 


10.909 (.02) 


10.763 


(.02) 


10.785 


(.02) 


837.2 


11.907 


(.04) 


10.968 (.04) 


10.794 


(.04) 


10.736 


(.04) 


839.2 


12.076 


(.06) 


11.009 (.06) 


10.795 


(.06) 


10.602 


(.06) 


844.3 


12.744 


(.07) 


11.432 (.07) 


11.113 


(.07) 


10.780 


(.07) 


849.2 


12.961 


(.02) 


11.762 (.02) 


11.395 


(.02) 


11.060 


(.02) 


853.2 


13.082 


(.06) 


11.972 (.06) 


11.670 


(.06) 


11.386 


(.06) 


856.6 


13.253 


(.21) 


12.085 (.21) 


11.805 


(.21) 


11.579 


(.21) 


862.7 


13.329 


(.02) 


12.242 (.02) 


12.052 


(.02) 


11.837 


(.02) 


867.2 


13.430 


(.02) 


12.452 (.02) 


12.230 


(.02) 


12.067 


(.02) 


871.6 


13.442 


(.02) 


12.505 (.02) 


12.341 


(.02) 


12.225 


(.02) 


872.7 


13.451 


(.02) 


12.571 (.02) 


12.374 


(.02) 


12.305 


(.02) 


877.6 


13.572 


(.10) 


12.714 (.10) 


12.537 


(.10) 


12.501 


(.10) 


881.7 


13.605 


(.08) 


12.840 (.08) 


12.661 


(.08) 


12.680 


(.08) 


888.6 


13.711 


(.12) 


12.932 (.12) 


12.904 


(.12) 


12.964 


(.12) 


889.6 


13.691 


(.05) 


13.032 (.05) 


12.932 


(.05) 


13.027 


(.05) 


894.7 


13.743 


(.05) 


13.145 (.05) 


13.095 


(.05) 


13.182 


(.05) 



Notes. JD - 2,455,000. Errors are given in parentheses. 



ing Cepheids jFreedman et all 1200 UlShappee & Stanekll201 Ih . 
Therefore, SN 201 Ife is an ideal object to test the current state- 
of-the-art of the SN la LC fitters. 

In this paper we present new, homogeneous, calibrated 
(BVRI) photometry for SN 2011fe obtained with a single tele- 
scope/detector combination (Sect. |2]l. We apply the two most 
widely accepted and trusted, ind ependently cahbrated, public 
LC fitters f or SNe la MLCS2k2 dJha. Riess & Kirshneit l2007h 
and SALT2 (iGuv et al.L 12007) to derive photometric distances to 
MlOl from our data (Sect. O. The results are compared with 
other MlOl distance estimates in Sect.|4] Section |5] summarizes 
our results. 



2. Observations 

2.1. Photometry of SN 2011 fe 

We have obtained multi-color ground-based photometric obser- 
vations for SN 201 Ife from the Piszkesteto Mountain Station 




E 

Fig. 1. The field around SN 201 Ife (marked by two line seg- 
ments) with the local comparison stars encircled. 

of the Konkoly Observatory, Hungary. We used the 60/90 cm 
Schmidt-telescope with the attached 4096 x 4096 CCD (FoV 
70x70 arcmin^, equipped with Bessel BVRI filters). In Table [1] 
the data for the first 41 nights are presented. Note that we use JD 
instead of MJD throughout this paper 

The magnitudes were obtained by applying aperture photom- 
etry using the daophot/phot task in /TJAiO Because the back- 
ground level around the SN position is relatively low and uni- 
form (see Fig.[T]i, neither PSF-photometry, nor image subtraction 
were necessary to get reliable light curves for SN 201 Ife. 

Transformation to the standard system was computed by us- 
ing color terms expressed in the following forms for the V mag- 
nitude and the color indices: 

V-v = Cv{V-I) + ^v 
(B-V) = Cbv ■ (b - v) + ^BV 
(V-R) = CvR ■ (v - r) + 

iV-I) = Cvi-(v-i) + (vi, (1) 

where lowercase symbols denote the instrumental magnitudes, 
while uppercase letters mean standard magnitudes. The color 
terms were determined by measuring Landolt standard stars in 
the field of PG2213 observed during photometric conditions: 
Cv = -0.019, Cbv = 1-218, Cvr = 1.035, Cvi = 0.959. These 
values were kept fixed while computing the standard magnitudes 
for the whole dataset. 

Zero-points (^x) for each night were measured using local 
tertiary standard stars (Table |2|. These local comparison stars 
were tied to the Landolt standards during the photometric cali- 
bration. 

Additional unfiltered photometry has been carried out at 
the Baja Observatory of Bacs-Kiskun County, Baja, Hungary 
with the 50 cm automated BART-telescope equipped with a 
4096 X 4096 back-illuminated Apogee Ultra CCD (FoV 40 x 40 
arcmin^, the frames were taken with 2x2 binning). During 

' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Table 2. Local tertiary standai-ds in the field of SN 2011fe 



ID 


RA (J2000) 


Dec (J2000) 


B (mag) 


V (mag) 


R (mag) 


I (mag) 


A 


14:04:04.438 


+54 


13 


32.64 


14.268 (0.028) 


13.404 (0.010) 


12.878 


(0.019) 


12.451 


(0.015) 


B 


14:03:45.175 


+54 


16 


16.31 


14.698 (0.033) 


14.054 (0.014) 


13.651 


(0.025) 


13.350 


(0.022) 


C 


14:03:28.982 


+54 


11 


33.74 


17.386 (0.127) 


16.320 (0.056) 


15.771 


(0.090) 


15.176 


(0.082) 


D 


14:03:24.941 


+54 


13 


57.36 


17.022 (0.117) 


16.405 (0.059) 


15.926 


(0.096) 


15.645 


(0.091) 


E 


14:03:23.779 


+54 


14 


32.83 


16.151 (0.068) 


15.617 (0.034) 


15.203 


(0.056) 


14.778 


(0.052) 


F 


14:03:22.410 


+54 


15 


36.22 


16.297 (0.055) 


14.910 (0.022) 


14.044 


(0.036) 


13.327 


(0.032) 


G 


14:02:38.490 


+54 


14 


50.69 


16.455 (0.085) 


16.012 (0.044) 


15.568 


(0.072) 


15.168 


(0.068) 


H 


14:03:05.865 


+54 


17 


25.49 


16.899 (0.101) 


16.160 (0.049) 


15.724 


(0.081) 


15.369 


(0.075) 


I 


14:03:05.803 


+54 


15 


19.91 


17.805 (0.162) 


16.625 (0.069) 


16.536 


(0.122) 


16.096 


(0.111) 


J 


14:02:54.159 


+54 


16 


29.17 


14.622 (0.032) 


14.043 (0.013) 


13.675 


(0.024) 


13.324 


(0.021) 



Notes. Magnitude errors are given in parentheses. 



Table 3. Unfiltered (scaled to /J-band) photometry of SN 201 Ife 
from the Baja Observatory, Hungary 



JD 


R (mag) 


J.D. 


R (mag) 


2455796.4 


> 18.50 (0.48) 


2455822.3 


10.21 (0.02) 


2455800.3 


13.11 (0.03) 


2455830.2 


10.77 (0.02) 


2455802.3 


12.03 (0.02) 


2455831.3 


10.81 (0.03) 


2455804.3 


11.25 (0.02) 


2455832.3 


10.81 (0.02) 


2455805.3 


10.96 (0.02) 


2455834.2 


10.88 (0.02) 


2455808.3 


10.41 (0.03) 


2455835.3 


10.89 (0.02) 


2455809.3 


10.27 (0.02) 


2455837.2 


10.95 (0.03) 


2455811.3 


10.10(0.02) 


2455838.3 


10.97 (0.02) 


2455815.3 


10.10(0.03) 


2455863.2 


12.16(0.02) 


2455817.3 


10.01 (0.02) 


2455866.2 


12.35 (0.03) 


2455818.3 


10.11 (0.02) 


2455867.2 


12.34 (0.03) 


2455819.3 


10.12(0.02) 


2455868.2 


12.32 (0.03) 


2455820.3 


10.16(0.03) 


2455871.7 


12.44 (0.02) 



Notes. Errors are given in parentheses. 
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Fig. 2. Comparison of LCs from Konkoly and RIT Observatories 



the course of the Baja-Szeged-Supernova-Survey (BASSUS) the 
field of MlOl was imaged with BART on 2011-08-22.9 UT, ~ 2 
days before discovery. No object was detected at the position of 
SN 20 life brighter than ~ 18.5 /?-band magnitude. After dis- 
covery, unfiltered photometric observations were taken on 25 
nights between Aug. 26 and Nov. 6, 2011 (TableO. These data 
were scaled to the properly calibrated 7?-band observations from 
Konkoly Observatory and used only in constraining the the mo- 
ment of explosion and the time of maximum light. 



Calibrated photometry for SN 201 Ife have al s o been col- 
lected from recent literature. [Richmond & Smithl (1201 2h pre- 
sented BVRI photometry obtained at the Rochester Institute 
of Technology (RIT) Observatory, and at the Michigan State 
University Campus Observatory. A comparison between the 
Konkoly and RIT data are plotted in Fig.|2](restricted to within 
40 days around peak for better visibility) illustrating the excel- 
lent agreement between these independent datasets. 

2.2. A UV-NIR spectrum of SN 2011 fe 

SN 201 Ife was intensively followed up spectroscopically, and 
the spe ctro scopic evolut i on is d iscussed in detail in Parrent et ^ 
(120 12h and ISmith et all (1201 lb . Expanding on that work is be- 
yond the scope of this paper Instead, we show only a single 
pre-maximum spectrum of SN 201 Ife (Fig. O extending from 
the ultraviolet (UV) to the near infrared (NIR). 

The spectrum plotted in Fig. [3] is a result of combining 
three datasets, obtained with different instruments. The opti- 
cal data were obtained with the HET Marcario Low Resolution 
Spectrograph (LRS, spectral coverage 4200 - 10200 A, resolv- 
ing power A/AA ~ 600) at the McDonald Observatory, Texas, 
on Aug. 27, 201 1. These data were reduced with standard IRAF 
routines. The UV part was taken b y Swifi/UVOJ as a U GRISM 
observation on Aug. 28, 2011 (see lBrown et al.Ll2012h . The re- 
duction was done with the routine uvotimgrism in HEASoft. 
The low {R ^ 200) and medium {R ^ 1200) resolution NIR 
spectra were obtained on August 26.3 UT with the 3.0 meter 
telescope at the NASA Infrared Telescope Facility (IRTF") usin g 
the SpeX medium-resolution spectrograph dRavner et al.Ll2003h . 
The IRTF data were reduced using a package of IDL routines 
spe cifically designed for the reduction of SpeX data (Spextool v. 
3.4. ICushing. Vacca. & Ravneiil2004t) . 

Fig. [3] illustrates the unprecedented quality of data available 
for SN 201 Ife, the analysis of which will be the topic of sub- 
sequent papers (e.g. a sequence of NIR spectra will be studied 
by Hsiao et al., in prep.). A similar extended spectr um for the 
Type la SN 201 liv has been recently published bv iFolev et alJ 
(l2012bh . The SN 201 Ife spectrum presented here is only the 
second such high-quality UVOIR spectrum for a Type la. These 
kind of data may be especially useful for theoretical modeling. 

Comparison with spectra of other SNe dCenko et al.L 1201 ih 
revealed that SN 201 Ife is a textbook-example of Branch- 
normal SNe la. The spectroscopic evidence that SN 201 Ife is 
a "normal" (i.e. not peculiar) SN la strengthens the applicability 
of the LC fitting techniques (see above) that were calibrated for 
such SNe. 
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SN' 2011 fe 2011-08-27 - 




3000 4000 5000 7000 10000 15000 20000 
Rest wavelength (A) 



Fig. 3. UV-NIR spectrum of SN 201 Ife. Note that scales on both 
axes are logarithmic. 

3. Analysis 

This section contains a brief review of the LC fitting methods for 
SNe la. Their application to the observed data of SN 201 Ife are 
then presented. 

3.1. SN la light curve fitters 

The empirical correlation between the LC shape and p eak bright- 
ness of SNe la was first revealed by Phillips ( il993h . after the 
initial suggestion made by Pskovskii C1977.) . According to the 
Phillips-relation, SNe la that dechne more slowly after maxi- 
mum are intrinsically brighter than the more rapidly declining 
ones. Phillips ( 1993) introduced the Aoti5(Z?) parameter for mea- 
suring the decline rate: it gives the decrease of the SN brightness 
from the peak magnitude at 15 days after maximum in the B- 
band. 

This concept was fu rther examined and extended to other 
photometric bands by iRiess et aL I (fT996l) . introducing the 
"Multi-Color Light Curve Shape" (MLCS) method. They defined 
a new parameter, A, for measuring the peak brightness as a 
function of the LC shape. Originally MLCS was calibrated for 
the Johnson-Cousins BVRI bands, and the LC in each band 
was described as a linear combination of two empirical (tabu- 
lated) curves and the parameter A. These curves were calibrated 
("trained") using 9 nearby, well-observed SNe la that had in- 
dependent distances, mostly from the Tully-Fisher method. In 
iRiess et alj (119981) MLCS was reformulated, expressing the LCs 
as a quadratic function o f A and including the U -band (but see 
e.g. iKessler etalll2009l for discussion on the utility of the U- 
band data). 

In this paper we appli ed the latest version, MLCS2k2 
(Qh a. Riess & KirshneR l2007h . which further improved the cali- 
bration by applying a sample of 133 SNe la for training and also 
included a new parametrization for taking into account the effect 
of interstellar reddening. In this version the observed LC of a SN 
la can be expressed as 

m,(f - to) = MO(f - to) +IUQ+ Ua, + §^)A^ + 

+ PAt - fo) ■ A + e,(f - to) ■ A\ (2) 

where t-to is the SN phase in days, to is the moment of maximum 
light in the B-band, m ^ is the observed magnitude in the x-band 
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Fig. 4. Fitting for the moment of explosion (left panel) and the 
time of Z?-band maximum (right panel) 

(x - B, V,R,I), M^l(t - to) is the fiducial la absolute LC in the 
same band, /zo is the true (reddening-free) SN distance modulus, 
f t, Cx and fSjc are functions describing the interstellar reddening, 
Ry and Ay are the ratio of total-to-selective absorption and the V- 
band extinction at maximum light, respectively, A is the main LC 
parameter, and P ^ and are tabulated functions of the SN phase 
("LC- vectors"). Together with A, the functions M*j*, P^ and 
describe the shape of the LC of a particular SN. For these func- 
tions we have applied the late st calibrat ion dow nloaded from the 
MLCS2k2 websiteB Note that lJha. Rie'ss & KirshneR (l2007l) tied 
the MLCS2k2 LC-vectors to SNe la in the Hubble-flow adopting 
Ho = 65 kms^'Mpc"'. Thus, if needed, the distances given by 
MLCS2k2 may be rescaled to other values of Ho by 

fioiHo) = A<o(MLCS)-51ogio(i?o/65)mag. (3) 

Another independent LC fitte r, SALT 2, was de veloped by the 
SuperNova Legacy Survey team ( Guv et ani2007h . SALT2 is dif- 
ferent from MLCS2k2 because it models the whole spectral en- 
ergy distribution (SED) of a SN la as 

FAp) = xo ■ [Mo(p, A) + xiMiip, A)] exp[c • CL(A)], (4) 

where p = f - fo is the time from Z?-maximum (the SN phase), 
is the phase-dependent rest-frame flux density, Moip, A), 
M\(p,A) and CL{A) are the SALT 2 trained vectors. The free pa- 
rameters Xo, x\ and c are the normalization- , stretch- and color 
parameters, respectively. We applied version 2.2.2b of the codqj, 
which was trained with the SNLS 3-year data (Guy et al., 2013, 
GIO hereafter). 

Because SALT2 models the entire SED, the observed LCs 
made with a particular filter set must be derived by synthetic 
photometry. SALT2 performs this computation based on the in- 
formation provided by the user on the magnitude system in 
which the input data were taken. Since our photometry is in the 
Johnson-Cousins system (see Sect.|2]i, we have selected the stan- 
dard Vega-magnitude system given in the code. 

3.2. Constraining the moment of explosion and B-band 
maximum 

Although the LC fitters applied in this paper use the time of B- 
band maximum light as the zero point of the time, the moment 

^ http : //www . physics . rutgers . edu/"saurabh/mlcs 2k2/] 
^ .http : //supernovae . in2p3 . fr/"guy/salt/usage . html] 
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of explosion is also a very important physical parameter for SNe 
la. This can be inferred from the pre-maximum photometry. In 
this section we repeat the analysis of Nugent et al. (20 11) to es- 
timate this parameter for SN 201 Ife using more data and a better 
sampled early LC. 

The pre-maximum LC of SNe la can be surprisingly well de- 
scribed by the constant temperature "fireball" m odel (e.g. lArnetH , 
Il982t iNugent et all 120111: iFolev et al.L l2012al) . In this simple 
model the adiabatic loss of the ejecta internal energy is just com- 
pensated by the energy input from the radioactive decay of the 
^*Ni and ^''Co synthesized during the explosion. This results in a 
nearly constant eff'ective temperature and a luminosity governed 
only by the change of the photospheric radius, which, assuming 
homologous expansion, can be approximated as Rph ~ t, giving 
L ~ fi. Although this simple picture does not capture all the de- 
tails in the pre-maximum eject a, it provides a surprisingly good 
fit to t he observations (see a lso Riess et anil999l;lHavden etall 
l2010HGaneshalingam et alll201 ih . 

We have applied the function of m = -2.5 logjo(A: ■ (f - texpf') 
to the observed pre-maximum /?-band magnitudes in Table [1] 
and|3] supplemented by similar data collected from Nuge nt et al.l 
(2011) and Richmond & Smith (2012). The fit parameters were 
k and texp, where the latter wa s further constrained by the epochs 
of pu blished non-detections dNugent et al.L 1201 ll: iBloom et al.L 
1201 2l and Table |3]l. The left panel of Fig. |4] shows the excel- 
lent agreement between the observed data and the fit law. 
The fitting resulted in f„p =JD 2,455,797.216 + 0.010, in per- 
fect agree ment with the value JD 2,455,797.187 + 0.014 re- 
ported by 'Nugen t et al.l (1201 ih . Relaxing the constraint to 
f and optimizing n gave n - 2.050 + 0.025 and texp - 
JD 2,455,797.182 ± 0.021, which do not diff^er significantly 
from the results assuming fi. We conclude that the explosion of 
SN 201 Ife occured at JD 2,455,797.20 + 0.16 (2011-08-23 
16:48 UT + 12 min). 

Th e B-band data from Table[T|and from [Richmond & SmithI 
( 1201 2|) are also used to constrain the moment of B-maximum. 
Fitting a 4th order polynomial to the magnitudes obtained be- 
tween -1-9 and -1-30 days after explosion resulted in tBmax - 
JD 2,455,814.4 + 0.6, which was used as input for the LC 
fitter codes (see below). Thus, the Z?-band maximum occured 
~ 17.2 days a f ter exp losion, very similar to the value derived 
bv lFolev et all ( 1201 2ab for SN 2009ig (17.13 days), which was 
also discovered in less than a day after explosion. The average 
value for the maior ity of "normal" SNe la is ~ 17.4 ± 0.2 days 
dHavden et all [201 Oi) . This supports the conclusion from spec- 
troscopy that SN 201 Ife is a "normal" SN la. 



3.3. Distance measurement 

We have applied both MLCS2k;2 and SALT2 to the observed BVRI 
data of SN 201 Ife shown in Sec.2.1. Before fitting, all data have 
been corrected for Milky Way reddening adopting Ay - 0.028 
mag and E(B - V) = 0.009 mag ( Schlegel et al.lll998 ). Note that 
these values are co nsistent with the recent recalibration of Milky 
Way reddening bv lSchlaflv & Finkbeine"3(l201 ih . 

Because of the low redshift of the host galaxy (z - 0.000804, 
see Sect. [T]i, A'-corrections for transforming the observed mag- 
nitudes to rest-frame bandpasses are negligible, th us, they were 
ignor ed. We have not included {/-band LC data (iBrown et al.L 
1201 2h in either fitting, thus avoiding the persistent systematic un- 
certa inties in modeling SNe la U -band data (e.g. iKessler et al.L 
12009ft . 
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Fig. 5. The fitting of the MLCS2k2 LCs to all observed data of 
SN 201 Ife. Solid curves represent the best-fitting templates, 
while dotted curves denote the template uncertainties given by 
the time-dependent variance of each template curve. 

Table 4. MLCS2k2 best-fitting parameters 



Parameter 


All data 


-7d < / < +40d 


to (JD) 


2,455,814.60 (0.10) 


2,455,814.9 (0.20) 


A (mag) 


-0.01 (0.08) 


0.02 (0.08) 


A^,"" (mag) 


0.05 (0.01) 


0.00 (0.02) 


yUo(//o = 73) (mag) 


29.21 (0.07) 


29.23 (0.07) 


Reduced x' 


0.2682 


0.2507 



Notes. Errors are given in parentheses. 



Note that the absolute magnitudes of SNe la LCs were cali- 
brated to different peak magnitudes in the two independent LC- 
fitters. MLCS2k2 was tied to SNe la dis t ances assuming Hq = 
65 km s 'Mpc~' dJha. Riess & Kirshned. l2007h . while the peak 
magnitude in SALT2 was fixed assuming Hq = 70 km s^'Mpc"' 
(GIO). To get rid of this discrepancy, we have transformed all 
distance moduli given below to Ho - 73 km s^'Mpc"' using 
Eq.m 

Since the MLCS2k2 templates are defined between -10 and 
+90 days around B-maximum, while the SALT2 templates ex- 
tend from -20 to only +50 days, we performed the LC fitting 
not only for all observed data (listed in Table [U, but also for 
those obtained between -7 and +40 days around B-maximum 
(JD 2,455,808 < t < 2,455,855). This test was performed in 
order to reach maximum compatibility between the applications 
of the two methods, and reduce the systematics that might bias 
the fitting results. 

3.3.1. MLCS2k2 

The fitting of Eq. |2] was performed by using a simple, self- 
developed -minimization code, which scans through the al- 
lowed parameter space with a given step and finds the lowest 
within this range. The fit parameters were the moment of B- 
maximum (fo), the V-band extinction Ay, the LC-parameter A 
and the distance modulus //q, with steps of 6to = 0.1, SAy = 0.01, 
dA = 0.01 and Sfio = 0.01, respectively. At the expense of longer 
computation time, this approach maps the entire hypersurface 
and finds the absolute minimum in the given parameter volume. 
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29.3 
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-0.1 0.1 
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Fig. 6. The map of the hypersurface around the minimum for 
MLCS2k2 fitting to all data. The contours correspond to 68, 90 
and 95 % confidence intervals (from inside to outside), respec- 
tively. 



We have fixed the reddening-law parameter as Ry = 3.1 
appropriate for Milky Way reddening, although several recent 
results suggest that some high-velocity S Ne la can be bet- 
ter modeled with s ignificantly lower Ry (IWang et al.L l2009t: 
lFolev&Kas"enLl20ilh . Since SN 201 Ife suff'ered from only mi- 
nor reddening and most of it is due to Milky Way dust (see be- 
low), it is more appropriate to adopt the galactic reddening law. 
Nevertheless, because of the low reddening, the value of Ry has 
negligible effect on the final distance. 

The best-fitting MLCS2k2 model LCs are plotted together 
with the data in Fig. |5] (we plot only the results of fitting the 
whole dataset, because the fit to the restricted data range gave 
very similar results). The final parameters are given in Table |4] 
for both the whole and the restricted data. The Icr uncertainties 
were estimated from the contour of A^^ - 1 corresponding to 
68% confidence interval. Fig.|6]shows the map of the the;^'^ hy- 
persurface and the shape of the contours around the minimum 
for the two key parameters A and jUq- It is seen that fio is strongly 
correlated with A, which is the major source of the relatively 
large uncertainty 6fi = 0.07 mag, despite the very good fitting 
quality. 

As seen in Table |4] there is no significant difference be- 
tween the fit parameters for the two datasets. The host extinction 
(A'y '^') is slightly less in the case of the restricted dataset, but 
that is compensated by the higher value of A (meaning fainter 
peak brightness), resulting in almost the same distance modulus. 
Thus, in the followings we adopt the parameters from fitting the 
full observed LC (left column in Table|4| as the final result from 
the particular LC-fitter, since those are based on the maximum 
available information. 

Note that the best-fitting iy[LCS2k 2 template LC corresponds 
to AmisiB) = 1.07 + 0.06. Although lRichmond & Sinithl(l2012h 
gives 1.21 + 0.03 for this parameter, this is probably a misprint 
since we measured 1.12 + 0.05 from t heir published data, simila r 
tothe value of 1.10 + 0.05 given bv Tam mann & Reindll (1201 ih . 
It seems that all these parameters are consistent with each other, 
as well as with the finding that SN 201 Ife has a nearly perfect 
fiducial SN la LC. 
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Fig. 7. The fitting of SALT2 LCs to 
SN 201 Ife. 
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3.3.2. SAU"2 

The SALT2 fitting was computed by running the code as de- 
scribed on the SALT2 website. The fit parameters provided by the 
code are: m* (rest-frame B-magnitude at maximum), xi (stretch) 
and c (color). Note that SALT2 restricts the fitting to data ob- 
tained no later than -1-40 days after maximum. 

Contrary to MLCS, the SALT2 model does not explicitly in- 
clude the distance or the distance modulus, thus, it must be de- 
rived from the fitting parameters. We followed two slightly dif- 
ferent procedures for this: th e one presented by G IO and an inde- 
pendent realization given bv Kessleretal] ( l2009i K09 hereafter). 
Starting from the fitting parameters m*, xi and c, the distance 
modulus fiQ in the GIO caUbration can be obtained as 

mBB = in*g - 0.008(+0.005) ■ xi + 0.013(+0.004) 
CC = 0.997(+0.097) ■ C + 0.002(±0.009) ■ xi + 0.035(±0.008) 

s = 0.107(+0.006)- jci -h0.991(+0.006) 
l-tQ - niBB - Mb + a ■ {s - V) - b ■ CC, (5) 

where we have adopted Mb - -19.218 ± 0.032, a - 
1.295 + 0.112 and/7 = 3.181 + 0.131 (GIO). 
The K09 calibration applies a simpler formula: 



/^o = - M() + a ■ x\ -/3 ■ c. 



(6) 



where Mq = -19.157 + 0.025, a = 0.121 + 0.027 and yS = 
2.63 + 0.22 have been adopted from K09. 

The uncertainties in the formulae above were taken into ac- 
count by a Monte-Carlo technique: we calculated 10,000 dif- 
ferent realizations of the above parameters by adding Gaussian 
random numbers having standard deviations equal to the uncer- 
tainties above to the mean values of all parameters and derived 
fio from each randomized set of parameters applying Eq. |5]and 
|6] Then the average and the standard deviation of the resulting 
sample of fio values are adopted as the SALT2 estimate for the 
distance modulus and its uncertainty. Table|5]lists the best-fitting 
parameters and errors, again for both the whole and the restricted 
dataset. The final SALT2 distance modulus was obtained as an 
unweigthed average of the two values from the GIO and K09 
calibrations. 
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Table 5. SALT2 best-fitting parameters 



Parameter 


All data 


-7d < / < +40d 


to (JD) 


2,455,815.505 (0.047) 


2,455,815.395 (0.097) 


ml (mag) 


9.959 (0.027) 


9.925 (0.029) 




-0.296 (0.049) 


-0.360 (0.080) 


c 


-0.030 (0.018) 


-0.060 (0.020) 


;Uo(G10) (mag) 


29.034 (0.078) 


29.088 (0.086) 


;Uo(K09) (mag) 


29.068 (0.062) 


29.105 (0.068) 


yuo(final) (mag) 


29.05 (0.08) 


29.10(0.08) 



Notes. Errors are given in parentheses. The final distance modulus (last 
row) is the average of the GIO and K09 estimates. 



4. Discussion 

The application of MLCS2k2 and SALT2 LC-fitters for all ob- 
served data resulted in distance moduli of /io(MLCS2k2) - 29.21 
+ 0.07 and /io(SALT2) = 29.05 + 0.08, respectively. It is seen that 
there is a ~ 2cr disagreement between these two values. Taking 
into account that these distance moduli were obtained by fitting 
the same, homogeneous, densely-sampled, high-quality photo- 
metric data of a nearby SN, and both codes provided excellent 
fitting quality, this A/io = 0.16 mag disagreement is rather dis- 
couraging. Note that the difference exists despite correcting the 
results from both codes to the same Hubble-constant, Hq - 13 
km s^'Mpc"' (see above). 

Restricting the LC fitting only to data taken between -7 and 
H-40 days around B-maximum (right column in Tables |4] and |5]l, 
the two distance moduli are both slightly higher and closer to one 
another: yL(o(MLCS2k2) = 29.23 + 0.07 and //o(SALT2) = 29.10 + 
0.08, giving AjUo = 0.13 mag. Since these parameters are gener- 
ally within the errors of those from fitting the complete LC, the 
~ 2cr disagreement still persists. 

A similar, even larger difference of 
//o(SALT2)-//o(MLCS2k2) = + 0.2 mag was found by K0 9 
for the "Nearby SNe" sample of iJha. Riess & KirshneH (l2007h . 
although deviations in both positive and negative directions 
have been revealed for individual SNe. Because the Nearby 
sample contains only a few very close, unreddened SNe like 
SN 201 Ife, the source of the mild discrepancy found by K09 is 
ambiguous. The present results suggest that, because of the lack 
of issues due to reddening, /T-coiTection, f/-band anomaly or 
spectral peculiarities in the case of SN 201 Ife (see Sect.[T]i, the 
> 0.1 mag difference between the MLCS2k2 and SALT2 distance 
moduli is probably entirely due to a systematic offset between 
the different zero-point calibrations of the fiducial SN peak 
magnitude in the two methods. 

In order to test this statement, we have compared the param- 
eters in Table|4]and|5] It is seen that for the whole dataset SALT2 
estimates the Z?-maximum (to) as being 0.9 day later than the 
epoch provided by MLCS2k2. The result from the simple polyno- 
mial fitting (Sect.lJSb is closer to the MLCS2k2 value, thus SALT2 
might tend to overestimate this parameter For the restricted data 
the final to from both methods changed slightly, becoming more 
similar, but a difference of ~ 0.5 day is still present. 

In order to investigate whether the uncertainty of to could be 
responsible for the systematic difference between the distance 
moduli, we have re-run the iy[LCS2k2 fitting for the restricted 
dataset by forcing to equal to the SALT2 value of 2,455,815.395. 
This resulted in A = 0.054 mag and fio - 29.205 mag with 
X'^ = 0.3242 (Ay'" remained the same). The changes of the pa- 
rameters are consistent with the shape of the surface plotted 
in Fig. |6l if A is increased, then /yq decreases. This lower fio is 



indeed closer to the SALT2 value, but a systematic difference of 
~ 0.1 mag still remains (i.e. the MLCS2k2 distance is still higher), 
while the quality of the fitting is clearly worse. Thus, while forc- 
ing f() to be equal to the SALT2 value may somewhat reduce the 
systematic difference between the two methods, not all of the 
systematics affecting the distance modulus can be explained by 
this parameter alone. 

Distance measurements independent from these SN la LC- 
fitters may help in resolving the open iss ue of absolute magni- 
tude and distance calibrations. Recently, Matheson et alJ (1201 2h 
published distance estimates to SN 201 Ife based on near- 
infrared (NIR) photometry. Because SNe la appear to be much 
better standard candles in the NIR than in optical bands, the 
usage of good-quality NIR photometry for this bright, nearby 
SN looks promising. Unfortunately, as iMatheson et alj (l20I2h 
concluded, the present state-of-the-art of getting SNe la NIR 
distances also suffers from an unsolved zero-point calibration 
problem. This resulted in a wide range of NIR distance moduli 
for SN 201 Ife spanning from 28.84 to 29.14 mag (corrected to 
Ho = 73, as above) with a mean value of yUo(NIR) - 29.0 + 0.2 
mag (Matheson et al., 2012). This is consistent with the SALT2 
result above, but it also agrees marginally (at ~ Icr) with the 
MLCS2k2 distance modulus. The range o f the NIR distance mod - 
uli from different calibrations, 0.3 1 mag (IMatheson et aUl20I2l) . 
is a factor of 2 larger than the uncertainties of the individual cal- 
ibrations (~ 0.15 mag) of NIR peak magnitudes of SNe la. This 
may also be a warning sign that the distance measurement tech- 
nique from NIR LCs of SNe la is far from being settled. 

The situation is not much better if one considers the vari- 
ous distance estimates available for the host galaxy, MlOl. This 
is one of the closest, brightest, and most thoroughly studied 
gala xies for which Cephei d-based distances are available (see 
e.g. IMatheson et al.L l20I2l and references therein). The most 
widely accepted distance modulus of MlOl is 29.13 + 0.11 
mag fro m the C epheid fL-relation by the HST Key Project 
( Freedman et aI.L 12001) . which is just in the middle between the 
MLCS2k2 and SALT2 distance m oduh above, being i n Icr ag ree- 
ment with both. More recently IShappee & Staned ( |20I I^ ob- 
tained 29.04 ± 0.19 mag from an independent study of MlOl 
Cepheids, which agrees better with the SALT2 estimate, al- 
though its larger errors makes the result also consistent with 
MLCS2k2. IShappee & Stanel3 (|20I Ih adopted the maser distance 
of NGC 4258 as their distance anchor, which gives a lower dis- 
tan ce modulus for t he LMC by 0.09 mag than the value adopted 
by iFreedman et al.l (t200lh . This accounts for most of the dif- 
ference between the two Cepheid-based results. Both of these 
Cepheid distances are slightly higher than the NIR-distance of 
Matheso n et al.l (l20I2h . but considering the larger errors of the 
latter (0.2 mag), the three estimates are all more-or-less consis- 
tent with each other 

Non-Cepheid distance estimates to MlOl span a ~ 0.35 
mag wide range, from 29 .05 (from the Tip of the Red Giant 
Branch method, TRGB, Shap pee & Stanekl l20IIh to 29.42 
(based on Planetary Nebulae Lumi nosity Function, PNLF, 
iFeldmeier. CiarduIIo & Jacobvl Il996h . which does not help 
m uch in resolving the is sue of the MlOl distance (see Fig. 3 
of IMatheson etalil20I2h . 

5. Conclusions 

The nearby, bright, weakly reddened Type la supernova 201 Ife 
in MlOl provides a unique opportunity to test both the precision 
and the accuracy of the extragalactic distances derived from SNe 
la LC fitters. In this paper we presented new, calibrated BVRI- 
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photometry for SN 201 Ife. The LCs were analyzed with publicly 
available LC-fitters MLCS2k2 and SALT2 to get the SN la-based 
distance to MlOl. There is a systematic offset of ~ 0.15 mag be- 
tween the MLCS2k2 and SALT2 distance moduli, the average of 
which also differs by ~ 0.13 mag from the distance estimate by 
iMatheson et"al] (1201 2h from SN 201 Ife NIR photometry. This 
systematic offset between the results of the two widely-used LC- 
fitters may be partly due to the different shape of the LC vectors 
near maximum, affecting the estimate of the moment of maxi- 
mum light, but the majority of the offset is probably caused by 
systematic errors of the peak magnitudes from different photo- 
metric calibrations. 

We conclude that the weighted average of the three dis- 
tance moduli of SN 201 Ife (using the inverse of the un- 
certainties as we ights), /zo(MLCS2k2), yL<o(SALT2) and jUo(NIR) 
dMatheson e t^Tl 2012h. and the two Cepheid-based dis tance 
to MlOl (iFreedman et al.L 120011; 'Shappee & Staneki 1201 lb pro- 
vides the following distance modulus of MlOl: 



/^O.MlOl 



29.109 + 0.049 (random) + 0.1 (syst) mag. 



(7) 



which corresponds to Dmioi = 6.6 + 0.5 Mpc, taking into ac- 
count both random and systematic uncertainties. Despite the ex- 
ceptional quality of the measured LCs of SN 201 Ife and the long 
history of efforts devoted to the calibration of LC-fitters for la 
SNe as well as the absolute distance to MlOl, the available tech- 
niques cannot predict the absolute distance to MlOl with better 
than 0.5 Mpc (~ 8%) accuracy. 
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